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25 

26 Abstract 

27 Background 

28 Understanding the evolution of the vertebrate pancreas is key to understanding its functions. The 

29 chondrichthyes (cartilaginous fish such as sharks and rays) have been suggested to possess the most 

30 ancient example of a distinct pancreas with both hormonal (endocrine) and digestive (exocrine) roles, 

31 although the lack of genetic, genomic and transcriptomic data for cartilaginous fish has hindered a more 

32 thorough understanding of the molecular-level functions of the chondrichthyan pancreas, particularly with 

33 respect to their "unusual" energy metabolism (where ketone bodies and amino acids are the main 

34 oxidative fuel source) and their paradoxical ability to both maintain stable blood glucose levels and 

35 tolerate extensive periods of hypoglycemia. In order to shed light on some of these processes we have 

36 carried out the first large-scale comparative transcriptomic survey of multiple cartilaginous fish tissues: 

37 the pancreas, brain and liver of the lesser spotted catshark, Scyliorhinus canicula. 

38 Results 

39 We generated a mutli-tissue assembly comprising 86,006 contigs, of which 44,794 were assigned to a 

40 particular tissue or combination of tissue based on mapping of sequencing reads. We have characterised 

41 transcripts encoding genes involved in insulin regulation, glucose sensing, transcriptional regulation, 

42 signaling and digestion, as well as many peptide hormone precursors and their receptors for the first time. 

43 Comparisons to published mammalian pancreas transcriptomes reveals that mechanisms of glucose 

44 sensing and insulin regulation used to establish and maintain a stable internal environment are conserved 

45 across jawed vertebrates and likely pre-date the vertebrate radiation. Conservation of pancreatic 

46 hormones and genes encoding digestive proteins support the single, early evolution of a distinct pancreatic 

47 gland with endocrine and exocrine functions in vertebrates, although the peptide diversity of the early 

48 vertebrate pancreas has been overestimated as a result of the use of cross-reacting antisera in earlier 

49 studies. A three hormone islet organ is therefore the basal vertebrate condition, later elaborated upon only 

50 in the tetrapod lineage. 
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51 Conclusions 

52 The cartilaginous fish are a great untapped resource for the reconstruction of patterns and processes of 

53 vertebrate evolution and new approaches such as those described in this paper will greatly facilitate their 

54 incorporation into the rank of "model organism". 
55 

56 Key words 

57 Catshark, pancreas transcriptome, Pdxl, Pdx2, insulin regulation 
58 

59 Background 

60 Chondrichthyans (cartilaginous fish such as sharks, skates, rays (elasmobranchs) and chimeras 

61 (holocephalans)) possess the earliest example of a distinct pancreatic gland containing multiple cell types 

62 with both endocrine and exocrine functions in vertebrates [1,2]. The more basal ("primitive") vertebrate 

63 lineages such as the jawless hagfish and lampreys (Figure 1) possess only small islet organs containing 

64 insulin- and somatostatin-producing endocrine cells and these islets lack any glucagon-producing cells or 

65 exocrine function [2-4]. The accumulation of multiple cell types into a single compact gland was an 

66 important step in pancreas evolution (and can be considered to be a vertebrate innovation) [5, 6] and it has 

67 been suggested that a switch from sensing gut-glucose to blood-glucose to establish a "stable inner milieu" 

68 via homeostatic mechanisms may have been an important factor in the evolution of a more complex 

69 glucose-dependent brain in vertebrates, protected from hyper- and hypoglycaemia [7-9]. However, the fact 

70 that insulin-like peptides in insects seem to fulfil similar roles in glucose metabolism and other 

71 physiological processes such as growth and reproduction suggests that at least elements of these 

72 mechanisms may have a more ancient origin [10]. 

73 Decades of research using light and electron microscopy and immunohistochemistry has revealed a great 

74 deal about the structure and organisation of the chondrichthyan pancreas and more recent studies have 

75 characterised the protein sequence and structure of some of the key pancreatic hormones [11-15]. The 

76 endocrine islets of chondrichthyans are typically scattered within exocrine tissue, sometimes associated 
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77 with minor ducts and contain a number of distinct cell types, thought to include the four major pancreas 

78 cell types: a-cells producing glucagon to increase blood glucose; (3-cells producing insulin to reduce blood 

79 glucose; 8-cells producing somatostatin to regulate pancreatic hormones such as insulin and glucagon and 

80 y-cells producing pancreatic polypeptide [16-18]. In addition to structural, cellular and hormonal 

81 conservation of the chondrichthyan pancreas compared to other vertebrates, there is also a conservation of 

82 function, with glucose-sensitive insulin release [11], pancreatectomy-induced hyperglycemia [19, 20] and 

83 exogenous insulin-induced hypoglycaemic effects [11, 21-23]. However, although blood plasma glucose 

84 levels are maintained at a fairly constant level during feeding and fasting (even over periods of up to 150 

85 days without food) [21, 24, 25], actual plasma glucose levels in elasmobranchs are lower than in teleost 

86 fish of comparative size and with similar metabolic rates [26]. It has also been found that elasmobranchs 

87 have an impressive tolerance of hypoglycaemia, including an ability to cope with a virtual absence of 

88 circulating glucose for at least 24 hours, a 75% reduction for at least a week and sub-normal plasma 

89 glucose levels for extended periods [22, 27] . There is an obvious paradox associated with an impressive 

90 ability to cope with long periods of hypoglycaemia existing in conjunction with the maintenance of 

91 apparently stable plasma glucose levels and others have pondered the necessity of central glucose-sensing 

92 mechanisms in these species [26]. 

93 The chondrichthyan pancreas represents an important model for studies of vertebrate pancreas evolution 

94 and function, particularly with reference to glucose homeostasis. However, full analysis of these areas has 

95 been hindered by a lack of genetic information and resources - almost the entirety of our current 

96 understanding of the chondrichthyan pancreas is based on what might be considered somewhat "old 

97 fashioned" (although still vital, important and informative) biological techniques, including descriptive 

98 gross anatomy and light and electron microscopy, enzymatic assays typically involving the injection of 

99 peptides derived from other (often mammalian) species and immunohistochemistry involving the use of 

100 antibodies raised against short mammalian peptide epitopes (see [2, 3, 26, 28] for reviews). There is 

101 currently a dearth of data regarding molecular level functions of the chondrichthyan pancreas, including 

102 mechanisms of transcriptional and translational control of gene regulation, signaling both in terms of cell- 
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103 cell communication within the pancreas and in terms of response to neuroendocrine signals and even more 

104 basic information such as the sequences of mRNA and protein precursors of previously identified or 

105 characterised digestive enzymes and pancreatic peptides. 

106 The lesser spotted catshark {Scyliorhinus canicula, often referred to as the lesser spotted dogfish) has 

107 recently become the chondrichthyan model of choice for a wide range of genetic, developmental and 

108 evolutionary analyses [29] and transcriptomic and genomic sequencing projects are currently underway 

109 for this species at Genoscope ( www.genoscope.cns.fr ). An enigmatic second member of the Pancreas and 

110 duodenal homeobox (Pdx) gene family (called Pdx2) was recently identified in the 5. canicula pancreas 

111 [30] but further studies of the role of this gene or the identification of the presence of additional members 

112 of other gene families involved in pancreas development, cell-specification and insulin regulation are 

113 impossible without more comprehensive molecular analyses. In order to shed further light on the possible 

114 role of the Pdx2 gene and to go some way to addressing the current dearth of data we set out to determine 

115 the pancreas transcriptome of the lesser spotted catshark and to carry out comparative expression analyses 

116 with other adult body tissues (liver and brain). These data represent the first large-scale transcriptomic 

117 analysis of multiple cartilaginous fish tissues and will be invaluable in understanding the functions of the 

118 cartilaginous fish pancreas, as well as shedding light on the evolution of the vertebrate pancreas itself. 
119 

120 Results 

121 A total of 6,260,398; 32,106,318 and 12,201,682 paired-end sequencing reads were generated for the 

122 pancreas, liver and brain respectively and these were pooled to generate a single assembly (Additional file 

123 1) which contained 86,006 contigs (when trimmed to remove all contigs <300bp, which likely represent 

124 single pairs of sequencing reads). The tissue distribution of these transcripts was determined by mapping 

125 sequencing reads from each tissue to this assembly and abundance values of >1 fragments per kilobase per 

126 million mapped reads (FPKM) were taken to confirm expression of a particular transcript in each tissue 

127 (Table 1). In this way 44,794 contigs were assigned to one or more tissues (Figure 2, Additional files 2-8). 

128 All transcripts contained an ORF encoding 20 amino acids or more (Figure 3), of which roughly 4-7% 
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129 encoded a signal peptide and so are likely to be secreted (Table 1). Typically, between 15-53% of 

130 transcripts had a BLAST hit in the RefSeq collection and 12-49% were annotated with GO terms (Table 

131 2). Although low, these figures are broadly comparable with a similar analysis of the white shark heart 

132 transcriptome [31], which found matches of 23.5% and 21.5% respectively. In order to provide a broad 

133 overview of the assigned gene ontology terms, we carried out a generic GOSlim annotation of the data 

134 (Figures 4 and 5), and Fishers exact tests showed that the pancreas was enriched for seven GO terms 

135 compared to liver ('cell', 'reproduction', 'transcription, DNA-dependent' , 'embryo development', 

136 'growth', 'intracellular non-membrane-bounded organelle' and 'non-membrane-bounded organelle') and 

137 two terms compared to brain ('transcription, DNA-dependent' and 'cellular amino acid metabolic process' 

138 - the full comparative enrichment results are provided in Additional file 9). 

139 A more detailed search strategy was carried out for particular categories of genes that would shed light on 

140 the similarities or differences of pancreas function in S. canicula compared to other vertebrates. The 

141 results of these analyses are outlined in the following sections. 

142 Pancreatic hormones and their receptors 

143 A large amount of immunohistochemical research has putatively identified the presence of a number of 

144 pancreatic peptide hormones in cartilaginous fish, including insulin, glucagon, somatostatin and pancreatic 

145 polypeptide, and the presence of at least some of these peptides has been confirmed by proteomic studies 

146 [11, 12, 15, 32]. Our transcriptomic data confirms the presence of mRNA transcripts encoding 

147 preproinsulin, preproglucagon (Figure 6) and preprosomatostatin (corresponding to the SSa gene [33]) in 

148 the pancreas and preprosomatostatin b and c in both the brain and pancreas, but not pancreatic polypeptide 

149 (PP - see next section), gastrin, gastric inhibitory polypeptide (GIP) or secretin. Vasoactive intestinal 

150 polypeptide (VIP) is expressed by both pancreas and brain, cholecystokinin by only brain and, as 

151 previously suggested [3], ghrelin appears to be absent from the shark pancreas and, indeed, from all three 

152 tissues. We find transcripts encoding the insulin receptor (IR) only in brain, the glucagon receptor only in 

153 liver and somatostatin receptor 1 (SSTR1) and SSTR5 in brain and both pancreas and brain respectively. 

154 Contrary to the findings of Larsson et al. [34], we find Neuropeptide Y receptors Yl, Y5 and Y6 in brain 
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155 only and Y8 in both pancreas and brain, suggesting that the expression of these receptors either varies 

156 between chondrichthyan species, or is more dynamic than previously assumed. 

157 The presence of PP and y-cells are key aspects of current schemes for the mode of vertebrate pancreas 

158 evolution [1-3, 35]. However, it has been known for some time that PP is tetrapod-specific, produced via 

159 duplication of the Peptide YY gene sometime prior to the divergence of this lineage [36, 37] and there is 

160 therefore a discrepancy between the findings of decades of immunohistochemical research and data from 

161 molecular genetic studies and analyses of vertebrate whole genome sequences. We have identified 

162 transcripts of two members of the Neuropeptide Y family (which includes Neuropeptide Y (NPY), Peptide 

163 YY (PYY) and Pancreatic polypeptide (PP)) in our dataset - a PYY gene expressed in pancreas and brain 

164 and a NPY gene expressed only in brain (Figure 7, both sequences are identical to published catshark 

165 sequences for PYY and NPY (accessions P69095 [38], AAB23237 [14]). We therefore suggest that older 

166 immunohistochemical studies which claimed to have detected PP+ cells in the cartilaginous fish pancreas 

167 may have in fact been relying on antisera that cross-reacted with PYY. A focus on the (often brief) 

168 methods sections of several key historical papers revealed that they in fact used the same anti-PP antibody, 

169 produced by Ronald Chance at Eli Lily in the 1970's [16, 39, 40]. It therefore appears that this antibody 

170 was detecting PYY in the pancreas of cartilaginous fish and that these initial papers and various 

171 subsequent papers have repeatedly been cited until the presence of PP in cartilaginous fish is considered to 

172 be established fact. In other cases, the misidentification of sequenced peptides has added to the confusion 

173 [38]. 

174 Our immunohistochemical surveys of the catshark pancreas using high-affinity anti-PP antibodies (Table 

175 1 in Additional file 10) showed varying results. The PP antisera from Sigma weakly stained the catshark 

176 pancreas but the staining was completely absorbed with PP, NPY and PYY peptides (Figure 8). While the 

177 Millipore anti-PP failed to stain (except on the mouse control pancreas, data not shown), experiments with 

178 anti-PYY antibodies detected strong signals, co-localising with insulin but not glucagon or somatostatin 

179 (Figure 8). NPY antisera immunoreactivity was detected in the same pattern as PYY (data not shown) and 

180 was absorbed with either PYY or NPY peptides (Table 2 in Additional file 10). 
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181 A search of the Blast2G0 results for the term 'Hormone activity' identifies three other peptide -encoding 

182 transcripts expressed in the catshark pancreas: Gastrin-releasing peptide (GRP, in pancreas and brain), 

183 which fulfils a variety of roles in the gastrointestinal tract, including the regulation of hormone release and 

184 the secretion of pancreatic enzymes [41, 42]; Neuromedin U (NMU, in all three tissues) which is 

185 expressed in nerves throughout the gastrointestinal tract [43] and Enkephalin (in pancreas and brain), an 

186 endogenous opioid that functions as a neurotransmitter or neuromodulator [44, 45]. 

187 Glucose sensing 

188 Hexokinase type IV, more typically called Glucokinase (GK), is a glucose-phosphorylating enzyme that 

189 has been shown to be the key molecule for glucose sensing in mammalian liver and pancreas cells [46] 

190 and mutations in GK are known to cause Maturity Onset Diabetes of the Young Type II (MODY2) [47, 

191 48]. Somewhat surprisingly, we find that GK is expressed in the shark brain and glucokinase regulatory 

192 protein (GKRP) only in liver. We also identified transcripts in all three tissues corresponding to 6- 

193 phosphofructo-2-kinase/fructose-2,6-bisphosphatase 2 {Pfkfbl), another regulator of glucose metabolism 

194 via interaction with GK. Finally we detected transcripts encoding two glucose transporters in all three 

195 tissues, namely the solute carrier family 2 genes Glutl and Glut2 and the shark pancreas therefore appears 

196 to represent an ancestral state where both types of transporter are expressed in the pancreas, as opposed to 

197 the the human pancreas (which relies on GLUT1) and that of rodents (which rely more heavily on Glut!) 

198 [49-51]. 

199 Insulin regulation 

200 The regulation of the insulin gene has, for obvious reasons, been an area of intensive study and a number 

201 of key regulators are now known from studies in rodents and humans. Perhaps one of the most important 

202 is the Pancreas and duodenal homeobox I (Pdxl) gene, also called Insulin promoter factor I (Ipfl), 

203 Islet/duodenum homeobox I (Idxl) or Somatostatin transactivating factor I (Stfl) [52-54] which (in 

204 addition to roles in embryonic development and (3-cell specification) binds to the TAAT motif-containing 

205 A boxes of the mammalian insulin promoter to stimulate transcription [55-58]. Mutations in PDX1 have 

206 been linked to Maturity Onset Diabetes of the Young Type IV (MODY4) [59] and pancreatic agenesis 
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207 [60, 61]. Cartilaginous fish and coelacanths have previously been shown to have retained an ancient 

208 paralogue of Pdxl, which we termed Pdx2 [30, 62] and accordingly we find transcripts of both genes in 

209 our pancreas transcriptome dataset. The basic helix-loop-helix (bHLH) transcription factor NeuroDl 

210 (previously called [32 [63]) has roles in pancreas development and islet formation [64, 65] and mutations 

211 in this gene have been linked to Type 2 Diabetes mellitus and MODY6 [66]. NeuroDl has also been 

212 shown to be important for insulin gene expression [67-69] and the LIM-homeodomain transcription factor 

213 Isll acts synergistically with NeuroDl and the bHLH transcription factor E47 to bind to and activate the 

214 insulin gene promoter [70-72]. The hepatocyte nuclear factor 1 alpha (Hnfla) gene was originally 

215 described as a liver-specific transcription actor, responsible for the regulation of a number of genes 

216 important for liver function [73-75]. However, it was later found that this gene also has a role in glucose 

217 homeostasis via regulation of insulin secretion [76, 77] and that mutations in Hnfla were the cause of 

218 MODY3 [78-80]. The Nkx6.1 homeodomain transcription factor is a potent transcriptional repressor with 

219 a key role in [3-cell differentiation [81, 82] and has also been shown to suppress the expression of 

220 glucagon to maintain [3-cell identity, as well as being able to regulate glucose-sensitive insulin secretion 

221 [83]. Our discovery of Pdxl (and Pdxl), NeuroDl and its partner E47, Isll, Hnfla and Nkx6.1 transcripts 

222 expressed in the catshark pancreas suggests an ancient role for these genes in vertebrate pancreas function 

223 and hints at early establishment of the insulin gene regulatory network. Additionally, the presence of 

224 transcripts encoding NeuroDl, e47, Isll and Nkx6.1 in the catshark brain highlights shared ancestry of 

225 these tissues in the vertebrate neuroendocrine system. We do not find any transcripts for MafA, which has 

226 been shown to be a key regulator of glucose-sensitive insulin secretion in humans and rodents [58, 84, 85], 

227 although other studies have also had difficulty identifying transcripts of this gene and other pancreas 

228 transcription factors in non-PCR based experiments [86, 87], possibly because of the low level of 

229 expression of transcription factors in general [88]. 

230 Transcription factors 

231 In addition to the transcription factors involved in insulin regulation discussed above, KEGG orthology 

232 (KO) analysis [89-91] has identified 13 transcription factors expressed in just the pancreas (including 
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233 Pdxl and Pdx2, FoxAl (Hnf3a) and Pancreas-specific transcription factor la (Ptfla)), 14 expressed in 

234 pancreas and liver, 33 in pancreas and brain and 51 in all three tissues (Tables 2 and 3). 

235 In a survey of the expression of 790 human DNA-binding transcription factors, Kong et al. [88] identified 

236 80 with expression restricted to the fetal pancreas, 32 restricted to the adult pancreas and 18 shared by 

237 both. Of the 31 adult-specific genes, we find evidence that 6 are also expressed in the adult catshark 

238 pancreas, although this number increases to 15 if members of the same gene family are considered (the 

239 possibility of divergent resolution of gene duplicates following the whole genome duplications [92] in 

240 early vertebrate ancestry must be considered). Since transcription factors are known to be expressed at low 

241 levels in cells (less than 20 copies per human adult cell [88]) it is likely that our figure is an underestimate 

242 and a more comprehensive survey of candidate transcription factor expression in this species is needed. 

243 Signalling 

244 Our KEGG orthology analysis identified 38 transcripts involved in signal transduction that are expressed 

245 only in the catshark pancreas, 1 1 in both pancreas and liver, 104 in pancreas and brain and 187 in all three 

246 tissues (Tables 2 and 3). Among these are representatives of the major vertebrate signalling pathways, 

247 including ligands and receptors for Fgf, Wnt, Notch, Vegf, Tgfp and Pdgf. Members of all of these 

248 pathways have previously been identified in the human pancreas transcriptome [87]. 

249 Homeobox gene diversity 

250 Homeobox genes are a group of transcription factors that encode a 60 amino acid DNA-binding 

251 homeodomain and that are involved in a wide variety of gene regulatory events in embryonic and adult 

252 tissues. A number of homeobox genes are known to be expressed during endodermal regionalisation and 

253 pancreas development, including Islet 1 and 2 (Isll, Isl2), Pancreatic and duodenal homeobox 1 {Pdxl), 

254 Nkx6.1, Nkx2.2, Pituitary homeobox 2 (Pitx2), Motor neuron and pancreas homeobox 1 (Mnxl), Onecut 

255 homeobox 1 (Onecut/Hnf6) and Paired box genes 4 and 6 (Pax4, Pax6) [93, 94]. Some older studies have 

256 detected a variety of homeobox genes in mammalian pancreas cell lines, including Cdx4, Hoxl.4 (HoxA4), 

257 Chox7 (Gbxl), Hox2.6 (HoxB4), Cdx3 (Cdx2), Cdxl, Hox4.3 (HoxD8), Hoxl.ll (HoxA2), Hox4a 

258 (HoxD3), Hoxl.3 (HoxA5) in the somatostatin-producing rat insulinoma cell line RIN1027-B2 [53] and 
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259 Isll, Lmx2, Alx3, HoxA4, HoxA13, Ipfl (Pdxl), Nkx2.2, Nkx6.1, En2 and Vdx in a hamster insulinoma cell 

260 line [95]. More recently, microarray and RNA-seq studies have identified a much larger number of 

261 homeobox genes expressed in the pancreas and especially the |3-ceil, with over 60 different homeobox 

262 genes identified by Kutlu et al. [87]. We used the homeodomain sequences of all human homeobox genes 

263 from HomeoDB [96, 97] and all vertebrate homeobox gene sequences from Pfam [98] as BLAST queries 

264 against our catshark transcriptome data and identified 1 1 different homeobox genes expressed in the 

265 pancreas, including five in just pancreas (HoxB5, HoxB7, Moxl, Pdxl, Pdxl), two in pancreas and liver 

266 (Hlx, Hhex), three in pancreas and brain (Arx, Zfhx3, Zfhx4) and one in all three tissues {Cut-like 2). These 

267 include genes known to be restricted to, or highly expressed in, [3-cells {Pdxl), a-cells {Arx) and acinar 

268 cell types {Cut-like 2) [86]. 
269 

270 Digestion 

271 In addition to its endocrine roles, the pancreas is also an important exocrine organ, fulfilling key functions 

272 in the digestion of proteins, lipids and carbohydrates. In the carnivorous elasmobranchs protein and lipids 

273 are the main energy sources [99] and it has been shown that ketone bodies and amino acids are the main 

274 oxidative fuel source for muscles and several other tissues, in preference to fatty acids [24, 28, 99]. 

275 Carbohydrates are thought to be utilised as oxidative fuels in elasmobranch heart muscle, as well as brain, 

276 red muscle and rectal gland [28, 100, 101]. It is therefore perhaps reasonable to assume that proteases and 

277 lipases are the most significant digestive enzymes produced by the elasmobranch pancreas and indeed this 

278 appears to be the case. Some form of chymotrypsinogen and trypsinogen have long been known to be 

279 produced by the elasmobranch pancreas, as has carboxypeptidase B, although these enzymes have not 

280 been fully characterised or isolated and sequenced [102-105]. We find transcripts of Elastase 2a and 3b, 

281 Chymotrypsinogen bl {Ctrbl), Chymotryp sin-like {Ctrl), Chymotryp sin-like elastase family, member 1 

282 {Celal) and Chymotrypsin-like elastase family, member 3B {Cela3b), Trypsin 1, 2 and 3, as well as the 

283 digestive carboxypepetidases (Al, A2, B 1) and those involved in activation and processing of other 

284 proteins, such as carboxypeptidase B2, D and E [106, 107]. 
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285 Some form of triacylglycerol lipase activity has previously been detected using crude enzyme preparations 

286 from the pancreas of skate (Raja (now Amblyraja) radiata) [108] and Leopard shark, Triakis semifasciata 

287 [109]. However, we find no evidence of Pancreatic liapse in the catshark pancreas transcriptome and 

288 instead find only Pancreatic lipase-related proteins 1 (Pnliprpl) and 2 (Pnliprp2). It is likely therefore that 

289 the triacylglycerol lipase activity found previously is a result of the action of Carboxyl ester lipase (CEL 

290 or bile salt-stimulated lipase) [1 10]. We also find colipase, in agreement with earlier studies of a range of 

291 cartilaginous fish and other basal vertebrates [111-113] and Hepatic and Hormone-sensitive lipases. 

292 Several lipid transporting apolipoproteins are also expressed by the catshark pancreas, including 

293 apolipoproteins A-IV, E, M and O. Finally, we have identified transcripts of genes involved in the 

294 digestion of carbohydrates (Pancreatic alpha amylase) and nucleic acids (deoxyribonuclease I and various 

295 ribonucleases). 
296 

297 Microsatellites 

298 It has recently been suggested [31] that a high frequency of dinucleotide simple sequence repeats (SSRs, 

299 microsatellites) is a general feature of shark genomes. We find 6,843 transcripts containing one or more 

300 di-, tri- or tetranucleotide microsatellites of five perfect repeats or more in our catshark data, with 482 of 

301 these only in pancreas, 3,083 only in brain and 473 only in liver (Table 4). In accordance with previous 

302 suggestions [31] we find dinucleotide repeats to be the most common type of SSR in both coding and non- 
303 coding regions of catshark transcripts. 

304 

305 Discussion 

306 Our analysis of the catshark pancreas transcriptome reveals the presence of genes known to be involved in 

307 glucose sensing and regulation of the insulin gene in other vertebrates and illustrates that functional 

308 conservation of these aspects of the vertebrate pancreas is reflected at the molecular -level. We therefore 

309 propose that these molecular-level mechanisms are a common feature of jawed vertebrates and that this 

310 lends support to the theory that the evolution of blood -glucose sensing and regulatory mechanisms may 
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311 have facilitated the evolution of the complex glucose-dependent brain of vertebrates [7-9]. We further 

312 suggest that the early evolution and fixation of these mechanisms has imposed evolutionary constraints on 

313 glucose sensing and insulin regulation in vertebrates, including in cartilaginous fish, even in the face of 

314 their ability to tolerate extended periods of hypoglycaemia and likely relaxed requirements for these 

315 processes. 

316 We find that the catshark pancreas produces at least eight peptide hormones (insulin; glucagon; 

317 somatostatin; peptide YY; gastrin-releasing peptide, neuromedin U, encephalin and vasoactive intestinal 

318 polypeptide, Table 5) and expresses a wide variety of genes involved in digestion, especially the digestion 

319 of proteins and lipids. The catshark pancreas therefore clearly has the features of a distinct pancreatic 

320 gland with both endocrine and exocrine functions and as such will be of great use in reconstructing the 

321 characteristics of the earliest vertebrate pancreas. The similarity in gene expression between the catshark 

322 and other vertebrates with respect to hormones, digestive enzymes, transcription factors and signaling 

323 pathways again provides support to the theory that there was a single, early origin of the pancreas at the 

324 base of the jawed vertebrate radiation. The overlap in peptides produced by the catshark pancreas and 

325 brain (Table 5) is a reflection of the shared ancestry of these tissues within the vertebrate neuroendocrine 

326 system [114]. 

327 Based on its co-localisation with insulin-, glucagon-, somatostatin- and PP-cells during mouse 

328 development, it has previously been suggested that a PYY+ cell may constitute a common progenitor of 

329 the major islet cell types [1 15]. Recent lineage tracing experiments have demonstrated that PYY+ cells 

330 give rise to islet d and PP cells and approximately 40% of pancreatic a and at least some P cells arose 

331 from peptide YY+ cells [116]. Most P cells and the majority of a cells are therefore not descendants of the 

332 peptide YY+/glucagon+/insulin+ cells that first appear during early pancreas ontogeny. The co- 

333 localisation of PYY with insulin in the adult shark pancreas illustrates the diversity of mechanisms that 

334 exist in vertebrate pancreas development and function and demonstrates the utility of "non-model" species 

335 to study these processes. The catshark PYY+ cells will therefore provide important insights into the 

336 evolution of the vertebrate pancreas, and especially progenitors of a, P, 8 and y-cells. 
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337 Our experiments make clear that much of the previous work on the presence or absence of peptides in 

338 basal vertebrate lineages may be suspect, with many false-positive signals resulting from cross-reacting 

339 antisera. Previous schemes of pancreas evolution based on these and similar data, which posited the 

340 restriction of various hormones to the alimentary canal (similar to the situation in protochordates such as 

341 amphioxus), the accumulation of these into a two-or three peptide islet organ in jawless fish and finally 

342 the "classic four-hormone islet tissue" of cartilaginous fish and other vertebrates [2] are therefore 

343 incorrect. In fact, it appears that the three hormone (glucagon, insulin, somatostatin) islet organ was 

344 established early in vertebrate evolution and remains today in the adult (but not larval) lamprey, 

345 cartilaginous fish and actinopterygian (ray-finned) fish, and that it is only in the sarcopterygian (lobe- 

346 finned fish) lineage that a four hormone (the above, plus PP) pancreas was formed. 

347 Our analysis of homeobox gene expression reveals a surprising level of variation between the genes 

348 known to be expressed in the catshark pancreas, human islets [87] and rat [53] and hamster [95] cell lines. 

349 It therefore seems likely that this particular class of transcription factors is extremely variable with respect 

350 to their spatial or temporal expression pattern in the vertebrate pancreas (or more likely both) and this is 

351 perhaps not too surprising given the variety of roles carried out by the pancreas in response to feeding, 

352 digestion and the regulation of blood glucose. As expected we have identified transcripts of both Pdxl and 

353 Pdx2 in the catshark pancreas, although we do not find any evidence for the presence of additional 

354 duplicates of other genes encoding proteins known to interact with PDX1 in other species. It therefore 

355 seems unlikely that the maintenance of paralogous Pdx2 genes in some vertebrate lineages reflects a wider 

356 conservation of duplicated gene regulatory networks produced as a result of whole genome duplication 

357 events early in vertebrate evolution. Comparison of the amino acid sequences of PDX1 and PDX2 across 

358 vertebrates shows conservation of the Pbx-interacting motif, DNA -binding domain and nuclear 

359 localisation signal but not of known transactivation domains and the PCIF1 -interaction domain [117-119] 

360 (Figure 9). The functions of the Pdx2 gene and the reasons for its retention in some species and 

361 independent loss in others (ray-finned fish and tetrapods) remain unknown. 



14 



Downloaded from http://biorxiv.org/on September 18, 2014 

362 With the availability of whole genome sequence information for a greater number of taxa and improved 

363 coverage of vertebrate pancreas transcriptomes a larger amount of data than ever before is now becoming 

364 available. These data, together with an appreciation that early vertebrate evolution was characterised by 

365 extensive genetic, developmental and morphological innovation facilitated by multiple whole genome 

366 duplications [120, 121] will better enable us to reconstruct pancreas evolution. As an example, we propose 

367 that the creation of the paralogous NPY and PYY during these duplications [34] facilitated the separation 

368 of the neuronal and gastroenteropancreatic (GEP) endocrine systems. We further suggest that the 

369 availability of additional copies of developmentally-important genes produced during the same duplication 

370 events [121] enabled the remodelling of the developing gut and the formation of a distinct pancreas with 

371 both endocrine and exocrine functions. 
372 

373 Conclusions 

374 We have generated a multi-tissue transcriptomic resource for an up and coming model organism, the 

375 lesser spotted catshark, Scyliorhinus canicula. Somewhat surprisingly we find few transcripts in common 

376 between the liver and pancreas, despite their relatively similar roles and shared developmental history as 

377 endodermal neighbors. The higher number of transcripts in common between brain and pancreas may 

378 provide evidence in support of the co-opting of neuronal programs by at least some pancreatic cells during 

379 vertebrate evolution [1 14, 122], although further comparative analyses are needed in this area. The 

380 similarity between the catshark pancreas transcriptome and those of various mammals with respect to 

381 insulin regulation, transcriptional and signaling machinery and peptide hormones and their receptors 

382 supports the single, early origin of a distinct pancreatic gland in vertebrates, although it seems likely that 

383 the peptide diversity of the early vertebrate pancreas may have been overestimated by older, 

384 immunohistochemical studies. The cartilaginous fish have a three peptide (insulin, glucagon and 

385 somatostatin) pancreas and the four peptide system seen in actinopterygian (ray-finned) and 

386 sarcopterygian (lobe-finned) fish and tetrapods is a later evolutionary innovation. The retention of the 

387 Pdx2 gene in cartilaginous fish does not apparently reflect a wider retention of duplicated members of 
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388 pancreas gene regulatory networks and the possible function(s) of this gene remains enigmatic. Our data, 

389 together with available or in progress transcriptomic and genomic resources for this and other 

390 chondrichthyan species will greatly facilitate comparative studies of elasmobranch, chondrichthyan and 

391 vertebrate evolution, particularly with reference to energy metabolism and the maintenance of stable blood 

392 glucose levels. 
393 

394 Methods 

395 RNA-Seq and sequence analysis 

396 Experimental methods involving animals followed institutional and national guidelines and were approved 

397 by the Bangor University Ethical Review Committee. Total RNA was extracted from freshly-dissected 

398 pancreas, liver and brain of two adult male and female catsharks approximately 24 hours post feeding. 

399 Pancreas samples were sequenced with using 2 x 250bp paired-end reads on the Alumina MiSeq platform 

400 at the Centre for Genomic Research (CGR) at the University of Liverpool. Brain and liver samples were 

401 sequenced using 2x 150bp paired-end reads on the Hlumina HiSeq 2000 platform at the Institute of 

402 Biological, Environmental & Rural Sciences (IBERS) at Aberystwyth University. Sequencing reads from 

403 the three tissues were assembled into a global tissue assembly using Trinity [123] with the jellyfish K-mer 

404 counting method. Tissue distribution of transcripts was assessed by mapping sequencing reads from each 

405 tissue to this global assembly, with an FPKM (fragments per kilobase per million mapped reads) value of 

406 >1 taken as confirmation of expression. Transcript annotation and assignment of gene ontology (GO) 

407 terms was performed using BLAST2GO [124, 125], the KEGG Automatic Annotation Server (KAAS 

408 [90]) and by local BLAST using BLAST+ v2.2.27 [126]. 
409 

410 Immunohistochemistry 

411 Male catsharks were euthanized according to a Schedule 1 method and the pancreas removed and fixed in 

412 4% paraformaldehyde/PBS overnight at 4°C. The fixed pancreas was then rinsed several times in PBS, 

413 dehydrated through a graded ethanol series and stored in 100% ethanol. 5pm sections of paraffin 
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414 embedded catshark pancreas were cut and mounted on glass slides. Slides were microwave treated in Tris- 

415 EGTA (TEG) buffer ph9.0 and allowed to cool for 30mins. Slides were rinsed in PBS and blocked in 

416 normal donkey serum and TNB buffer (Perkin Elmer) for 30mins. Primary antisera were added over night 

417 at room temperature and the next day the slides were rinsed 3 x 5min each in PBS and specific cross 

418 absorbed donkey anti- mouse, rabbit, or guinea pig secondary antisera (Jackson Immunoresearch) were 

419 added for 30mins. The slides were rinsed in PBS and mounted. Details of antisera are given in Table 1 in 

420 Additional file 10. All pictures were taken on a Zeiss Meta510 confocal microscope. 
421 

422 Antibody Absorption 

423 In order to test their specificity against the Pancreatic polypeptide family, the antisera were incubated 

424 overnight at 4°C with lOug of either pancreatic polypeptide (Sigma), Neuropeptide Y (Bachem) or 

425 peptide YY (in-house synthesis) or no peptide. The next day the antisera were added to the slides and the 

426 staining was performed as above. The staining intensity was compared to the no peptide control and given 

427 a rating of 1-3 (+, ++, +++). The results are shown in Table 2 in Additional file 10. 
428 
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774 Figure legends 

775 

776 Figure 1. Phylogenetic tree of the major extant vertebrate groups. The relationships of the most common 

777 chondrichthyan (cartilaginous fish) model species (Elephant shark, Callorhinchus milii; Little skate, 

778 Leucoraja erinacea; Lesser spotted catshark, Scyliorhinus canicula; Spiny dogfish, Squalus acanthias) are 

779 shown, as are representative lineages from the ray-finned (actinopterygian) and lobe-finned 

780 (sarcopterygian) fish. The origin of the combined endocrine and exocrine pancreatic gland at the base of 

781 the jawed vertebrates is indicated. 
782 

783 Figure 2. Tissue distribution of transcripts, as determined by mapping the sequencing reads derived from 

784 each tissue to a combined, all -tissue assembly. Contig values of >1 FPKM (Fragments Per Kilobase of 

785 exon per Million fragments mapped) were taken as evidence for expression. 
786 

787 Figure 3. Length and tissue distribution of open reading frames (ORFs) derived from assembled contigs 
788 

789 Figure 4. Proportion of transcripts assigned to each of the top 25 gene ontology (GO) slim 'Biological 

790 Process' terms for catshark pancreas, brain and liver tissue-specific transcripts. 
791 

792 Figure 5. Proportion of transcripts assigned to each of the top 25 gene ontology (GO) slim 'Molecular 

793 Function' terms for catshark pancreas, brain and liver tissue-specific transcripts. 
794 

795 Figure 6. Annotation of the precursor peptides of catshark preproinsulin and preproglucagon. Signal 

796 peptides (amino acids 1-24 and 1-20 respectively) are underlined and basic amino acid cleavage sites are 

797 lowercase. Glucagon -like peptides (GLP) la and lb are annotated based on similarity to the duplicated 

798 GLP1 peptides in the unpublished Squalus acanthias and Hydrolagus colliei proglucagon sequences 
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799 available on Genbank (accession numbers AAS57653 and AAS57654). An oxyntomodulin-like peptide 

800 has been purified from H. colliei and corresponds to amino acids 47-82 in the Catshark preproglucagon. 
801 

802 Figure 7. Amino acid alignment of vertebrate Peptide YY (PYY), Neuropeptide Y (NPY) and Pancreatic 

803 polypeptide (PP) sequences. Genbank accession numbers are given in square brackets. Sea, Scyliorhinus 

804 canicula (lesser spotted catshark); Sac, Squalus acanthias (spiny dogfish); Ler, Leucoraja erinacea (little 

805 skate); Cmi, Callorhinchus milii (elephant shark); Hsa, Homo sapiens (human); Lfl, Lampetra planeri 

806 (brook lamprey); Loc, Leucoraja ocellata (winter skate). 
807 

808 Figure 8. Immunolocalization of pancreatic hormones and pancreatic polypeptide and PYY in catshark 

809 pancreas. (A) The distribution of the pancreatic hormones insulin (blue), glucagon (green) and 

810 somatostatin (Red) in uniquely shaped islet structures. (B) Pancreatic polypeptide (red) specific antisera 

811 fail to stain a specific subset of endocrine cells in the pancreas, while insulin (blue) and glucagon show a 

812 normal distribution. (C-D) PYY shows colocalization with most of the insulin immunoreactive cells but 

813 not glucagon or somatostatin. All images are 250x magnification 
814 

815 Figure 9. Protein domains in vertebrate PDX1 and PDX2. Transactivation domains A-E [1 19], PCIF1- 

816 interaction domains [118], homeodomains, DNA-binding domains (i) and nuclear localisation signals (ii) 

817 are highlighted. Domain E contains the PBX-interacting hexapeptide motif [117]. There is very little 

818 conservation of amino acid sequence between the paralogous PDX1 and PDX2 suggesting that they carry 

819 out distinct functions within the pancreas, although clearly both are localised to the nucleus, bind DNA 

820 and interact with PBX proteins. Hsa, human {Homo sapiens); Mmu, mouse (Mus musculus); Rno, rat 

821 (Rattus norvegicus); Gga, chicken {Gallus gallus); Acar, Anole lizard (Anolis carolinensis); Xla, Xenopus 

822 laevis; Xtr, Xenopus tropicalis; Lme, Indonesian coelacanth (Latimeria menadoensis); Lch, African 

823 coelacanth {Latimeria chalumnae); Acal, Bowfin (Amia calva); Loc, Spotted gar (Lepisosteus oculatus); 

824 Dre, zebrafish (Danio rerio); Tru, fugu (Takifugu rubripes); Ola, medaka (Oryzias latipes); Gac, 
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825 stickleback; Tni, Green spotted puffer (Tetraodon nigroviridis); Sea, lesser spotted catshark (Scyliorhinus 

826 canicula); Ler, little skate (Leucoraja erinacea); Cmi, elephant shark {Callorhinchus milii); Bfl, 

827 amphioxus (Branchiostomafloridae) 
828 
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863 Table 1. Assembly characteristics. Tissue distribution of transcripts was assigned by read mapping, taking 

864 a value of >1 fragments per kilobase per million mapped reads (FPKM) as evidence of expression. 



Tissue(s) 


Number of 

/■•/■VfiTi rvc 

COIlLlgS 


Contig 

(hr>\ 
\ u t>) 


Longest 

COIlllg 

fhnl 


Mean ORF 
length (bp) 


IViax L/ivr 

leiigui 

fhnl 
\ u t>) 


Willi Slgllal 


L o.lll-'lCa.a Ullly 




49Q 


JJOU 


177 

1 / / 




9S7 n Q9%"i 


Rt*ain rvnlv 
ULtxm xJiiLy 




7Q8 


1 (S085 


338 

JJO 


1 3751 


1f)9Q C4 84%1 


T ,\ vpr onl v 


3423 


676 


13906 


271 


13532 


220 (6 43 %) 


Pancreas, brain 


3847 


1479 


10710 


482 


10394 


210 (5.46%) 


Pancreas, liver 


924 


1039 


5907 


329 


4505 


62 (6.71%) 


Brain, liver 


2482 


1329 


12080 


481 


10403 


149 (6.00%) 


Pancreas, liver, brain 


6326 


2257 


11135 


705 


9329 


356 (5.63%) 



865 
866 



43 



Downloaded from http://biorxiv.org/on September 18, 2014 



867 Table 2. Number of contigs with a BLAST hit in the RefSeq database and Gene Ontology (GO) 

868 annotation assigned by BLAST2GO [124, 125]. The number of transcription factors and signaling 

869 pathway components in each tissue or tissue combination is also shown. 
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874 Table 3. Transcription factors and signaling pathway components unique to the lesser spotted catshark 

875 pancreas based on available data. 



Pancreas-specific transcription factors 
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Nuclear receptor subfamily 0 group B 1 


Epidermal growth factor receptor 


Pancreas-specific transcription factor la 


Fibroblast growth factor 


Transcriptional enhancer factor 


FMS-like tyrosine kinase 1 


Zinc finger protein GLI3 


Frizzled 9/10 




Glutamine synthetase 




Inositol 1,4,5-triphosphate receptor type 1 




Inositol 1,4,5-triphosphate receptor type 2 




Inositol 1,4,5-triphosphate receptor type 3 




Insulin 




Integrin alpha 1 




Integrin alpha 2 




Integrin beta 6 
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Inter leukin 10 receptor beta 




Janus kinase 2 




Laminin, alpha 3/5 




Laminin, beta 4 




Leucine-rich repeats and death domain-containing protein 




Mitogen-activated protein kinase kinase kinase kinase 4 




Nucleoprotein TPR 




Protein crumbs 




Receptor-interacting serine/threonine-protein kinase 1 




Secretory phospholipase A2 




Suppressor of cytokine signaling 1 




TGF-beta receptor type-2 




Transcriptional enhancer factor 




Transferrin 




Vascular endothelial growth factor C/D 




Zinc finger protein GLI3 
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882 Table 4. Predicted di-, tri- and tetranucleotide simple sequence repeats (microsatellites) in each catshark 

883 tissue or combination or tissues. Results are shown for both full transcripts and predicted open reading 

884 frames (ORFs) and in both cases dinucleotide repeats are the most common. 





Transcript 


Predicted ORFs 




Di- 


Tri- 


Tetra- 


Di- 


Tri- 


Tetra- 


Pancreas only 


456 
(6.96%) 


23 
(0.35%) 


3 

(0.05%) 


60 
(0.92%) 


6 

(0.09%) 


0 

(0%) 


Brain only 


2808 
(13.22%) 


268 
(1.26%) 


7 

(0.03%) 


829 
(3.90%) 


147 
(0.69%) 


1 

(<0.01%) 


Liver only 


446 
(13.03%) 


24 
(0.70%) 


3 

(0.09%) 


168 
(4.90%) 


12 

(0.35%) 


1 

(0.03%) 


Pancreas + 
Brain 


481 
(12.50%) 


43 
(1.12%) 


2 

(0.05%) 


72 
(1.87%) 


15 

(0.39%) 


2 

(0.05%) 


Pancreas + 
Liver 


105 
(11.36%) 


5 

(0.54%) 


1 

(0.11%) 


44 

(4.76%) 


3 

(0.32%) 


1 

(0.11%) 


Brain + Liver 


702 
(28.28%) 


61 

(2.46%) 


3 

(0.12%) 


106 
(11.47%) 


21 

(2.27%) 


0 

(0%) 


Pancreas + 
Liver + Brain 


1240 
(19.60%) 


158 
(2.50) 


4 

(0.06%) 


308 
(4.87%) 


108 
(1.71%) 


2 

(0.03%) 
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890 Table 5. Peptide diversity of the catshark pancreas, brain and liver. Our comprehensive transcriptomic 

891 survey of the lesser spotted catshark pancreas highlights the disparity in the estimation of peptide diversity 

892 in early vertebrates as previously suggested by immunohistochemical (IHC) studies and highlights the 

893 similarity of pancreas and brain peptide complements. 



Peptide 


Pancreas 
IHC 


Pancreas 
transcriptome 


Brain 
transcriptome 


Liver 
transcriptome 


Insulin 


+ 


+ 


- 


- 


Glucagon 


+ 


+ 


- 


- 


Somatostatin 


+ 


+ 


+ 


- 


Pancreatic polypeptide 


+ 


- 


- 


- 


Peptide YY 


+ 


+ 


+ 


- 


Gastrin-releasing peptide 


+ 


+ 


+ 


- 


Neuromedin U 


+ 


+ 


+ 


+ 


Enkephalin 


+ 


+ 


+ 




Cholecystokinin 


+ 




+ 




Gastrin 


+ 








Vasoactive intestinal polypeptide 


+ 


+ 


+ 




Gastric inhibitory polypeptide 


+ 








Secretin 


+ 
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898 Additional files 

899 Additional file 1 : Combined tissue assembly, trimmed to remove contigs <300bp 

900 Additional file 2: Pancreas-specific transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

901 Additional file 3: Brain-specific transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

902 Additional file 4: Liver-specific transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

903 Additional file 5: Pancreas/brain transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

904 Additional file 6: Pancreas/liver transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

905 Additional file 7: Brain/liver transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

906 Additional file 8: Pancreas/brain/liver transcripts present at >1 FPKM, trimmed to remove contigs <300bp 

907 Additional file 9: Gene ontology (GO) enrichment results for pairwise tissue comparisons 

908 Additional file 10: Antibody table and peptide absorption results 
909 

910 
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911 Additional fde 10. 

912 Additional table 1. Table of antisera used in immunohistochemical surveys of the catshark pancreas. 

913 PYY, peptide YY; NPY, neuropeptide Y: PP, pancreatic polypeptide. 



Antigen 


Species 


Dilution 


Source 


PYY 


Rabbit 


1:500 


Sigma 


PYY 


Rabbit 


1:500 


AbCam, UK 


NPY 


Rabbit 


1:1000 


AbCam,UK 


Insulin 


Mouse 


1:75 


Novo Nordisk A/S, Denmark 


Insulin 


Guinea Pig 


1:200 


DAKO, Glostrup, Denmark 


Glucagon 


Mouse 


1:50 


Novo Nordisk A/S, Denmark 


Glucagon 


Rabbit 


1:500 


DAKO, Glostrup, Denmark 


PP 


Rabbit 


1:500 


Sigma 


PP 


Guinea Pig 


1:500 


Linco/Millipore, , 


Somatostain 


Mouse 


1:100 


Novo Nordisk A/S, Denmark 
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922 Additional table 2. Absorption of Pancreatic Polypeptide Family antisera. Staining is characterised from 

923 strong (+++) to weak (+) or absent (-). PYY, peptide YY; NPY, neuropeptide Y: PP, pancreatic 

924 polypeptide. 



Antibody 


Peptide 


Staining Dogfish 


Staining Mouse 


PP 


None 


+++ 


- 


PP 


NPY 


- 


- 


PP 


PYY 


- 


- 


PP 


PP 


- 


- 


PYY abeam 


none 


++ 


++ 


PYY abeam 


PP 


++ 


++ 


PYY abeam 


NPY 


+ 


++ 


PPY abeam 


PYY 


- 


- 


PYY Sigma 


none 


+++ 


+ 


PYY Sigma 


PP 


+++ 


+ 


PYY Sigma 


NPY 


+++ 


+ 


PYY sigma 


PYY 






NPY 


none 


+++ 


+++ 


NPY 


PP 


+++ 


+++ 


NPY 


PYY 






NPY 


NPY 
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